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1. Introduction 

The concept of quantum mechanics on a graph is more than half a century old 
having roots in modelling of aromatic hydrocarbons pQ. For many years, however, 
it was rather a curiosity, or maybe an interesting textbook example. The situation 
changed two decades ago with the advent of microfabrication techniques which allow 
us to produce tiny graph-shaped structures of semiconductor and other materials 
for which this is a useful and versatile model. This motivated a new theoretical 
attention to the subject - see, e.g., |2|3j . Since then the literature on quantum 
graphs grew to a formidable volume, and we restrict ourselves here to mentioning 
recent reviews in [4 5j5] where an extensive bibliography can be found. 

From the mathematical point of view the attractive feature of the model is that 
it deals with families of ordinary differential equations, the solutions of which have 
to be properly matched at the graph edge endpoints. Since the solutions are often 
explicitly known, the spectral analysis can be reduced to an algebraic problem. 
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The key point here are the boundary conditions through which the wave func- 
tions are matched. The Hamiltonian is typically a second-order differential operator, 
for instance, in the simplest case of a free spinless particle it acts on the j-th. edge as 
Hipj — —ipj. Thus the boundary conditions are linear relations coupling the values 
of the functions and their first derivatives at graph vertices; from the physical point 
of view it is usually sufficient to consider only local couplings which involve values 
at a single vertex only. Another general physical restriction is the self-adjointness of 
the Hamiltonian; it implies that a vertex joining n graph edges may be characterized 
by boundary conditions involving n 2 real parameters [3 3 . 

This leaves a considerable freedom in the choice of a model to describe particular 
physical systems, and an understanding of the physical meaning of vertex coupling is 
needed to pick the appropriate operator from the class of admissible Hamiltonians. 
A natural way to approach this problem is through approximation, i.e. regarding 
the quantum graph in question as a limit of a family of more "realistic" systems 
with a less number of free parameters. One possibility is to approximate a graph 
by a family of "fat graphs" or similar manifolds equipped with the corresponding 
Laplace-Beltrami operators. The best studied case is the one where the approxi- 
mated manifolds have Neumann boundary, or no boundary et all I? 8 I9I10I11I12) . 
where unfortunately the limit yields - of the multitude of available boundary con- 
ditions - only the most simple ones. There are also fresh results [13114) on the case 
with Dirichlet boundary but in general the approach based on squeezed manifolds 
did not yield so far a satisfactory answer to the question. 

Another, maybe less ambitious approach is to model vertex boundary conditions 
through families of interactions on the graph itself. Here two cases have to be 
distinguished. In the n 2 -parametcr family mentioned above the boundary conditions 
with wave functions continuous at the vertex form just one-parameter subfamily. 
These boundary conditions can be approximated by families of scaled potentials in 
analogy is analogy with one-dimensional S interactions [15] . The remaining, more 
singular cases require a different approach. An inspiration may be derived from the 
approximation of one-dimensional 8' interactions suggested, somewhat surprisingly, 
by Cheon and Shigehara in [16] and elaborated in a mathematically consistent way 
in |17|18j . It is based on a family of 8 interactions which approach each other being 
scaled in a particular nonlinear way. An analogous procedure for vertices of degree 
n > 2 was proposed in [19] in the case of the so-called S' s coupling; the key element 
here was the symmetry with respect to permutation of the edges which allowed 
to reduce the analysis to a one-dimensional halfline problem. The same technique 
was afterwards in [20] applied to the class of all permutation-symmetric boundary 
conditions which form a two-parameter subfamily in the n 2 -parameter set. 

The main goal of the present paper is to explore whether the idea of [16] can 
be adapted to situations without a permutation symmetry and how wide class of 
boundary conditions can be in this way described. As in the work mentioned above 
we will consider a star graph with a single vertex and n semi-infinite edges. For 
simplicity we will also assume that the motion on graphs edges is free; the obtained 
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approximations extend easily to Schrodinger operators on the graph provided the 
potentials involved are sufficiently regular around the vertex. We are going to show 
that the Cheon-Shigehara technique can produce for n > 2 at most a 2n-parameter 
family of boundary conditions at the vertex. Furthermore, we will demonstrate that 
such an approximations, with two 5 interaction at each edge, do indeed exist and 
that they converge in the norm resolvent topology. 

The next question is how to extend the approximation to a wider class of cou- 
plings. A natural possibility is amend the star by extra edges supporting S interac- 
tions which shrink to the "main" vertex with the parameter controlling the approx- 
imation. We devise such a scheme a show that it yields an ("J 1 ) -parameter family, 
generically all couplings which are time-reversal invariant. In this case, however, 
we restrict ourselves to deriving the boundary condition formally. We are convinced 
that the norm resolvent convergence could be verified as in the case mentioned but 
the argument would be extremely cumbersome. Notice that the idea of using ad- 
ditional edges to model singular couplings appeared already in 21 1. In contrast to 
that paper, however, we keep here the number of added edges fixed. 

Let us review briefly the contents of the paper. In the next section we gather 
the needed preliminary information. We review the quantum graph concept, recall 
different vertex couplings and review briefly the known approximations. In Section 3 
we analyze a CS-type approximation to the vertex in a star graph based on adding 
8 interactions on star edges, the following section is devoted to the proof of norm- 
resolvent convergence. Finally, in Section 5 we will describe the mentioned more 
general approximation with extra edges added to the star graph. 

2. Preliminaries 
2.1. Quantum graphs 

Let us first recall a few basic notions. A graph T is an ordered pair V = (V, E), 
where V and E are finite or countably infinite sets of vertices and edges, respectively. 
Without loss of generality we may identify E with a family of two-element subsets 
in V, excluding thus loops and multiple edges, since in the opposite case we can 
simply add extra vertices. The vertex degree of v G V is the number of edges which 
have v as its endpoint. T is a metric graph if each of its edges can be equipped with 
a distance, i.e. identified with a finite or semi-infinite interval of length I G (0, +oo]; 
the endpoints "at infinity" are conventionally not counted as vertices. In particular 
a star graph has a finite number n > 2 of edges and a single centre which is the 
only vertex where all the edges (called also arms in this case) meet. 

The subject of our interest is quantum mechanics on graphs. Given a metric 
graph r with edges J\, . . . , J n we identify the orthogonal sum H = ©™ =1 L 2 (Jj) 
with the state Hilbert space, i.e. the wave function of a spinless particle "living" 
on r can be written as the column \I/ = (ipi,ip2, ■ ■ ■ ,ipn) T with ipj G L 2 (Jj). In 
the absence of external fields the Hamiltonian H acts as (Hr x &)j = —ip", where 
as usual we put h = 2m = 1. Its domain consists of functions from W 2,2 (T) :— 



7, 2008 15:39 WSPC/INSTRUCTION FILE vertapp 



4 P. Exner and O. Turek 

d)J=i ^ 2 ' 2 (Jj)i since H is required to be a self-adjoint operator they must satisfy 
appropriate boundary conditions at the vertices which we will recall below. 

The meaning of these boundary condition is our main concern in this paper, 
therefore we restrict ourselves to graphs with a single vertex, namely star graphs 
with n semi-infinite edges Jj ~ R + , j = 1, . . . , n; we denote them as r or F„. 

2.2. Vertex couplings 

Since the Hamiltonian mentioned above is a second-order operator, the matching 
conditions involve boundary values of the functions in the vertex and of their first 
derivatives. Both regarded as one-sided limits, the derivatives are taken in the 
outward direction. We arrange them into column vectors ^(0) and ^'(0). The self- 
adjointness of H, which in the physical language means conservation of probability 
current at the vertex, is expressed through a linear relation between these vectors, 

A*(0) + £#'(0) = , (2.1) 

by [22] the operator H is self-adjoint if and only if A, B € C n,n satisfy the conditions 

rank(A, B) = n , AB* is self-adjoint, (2.2) 

where (A, B) denotes the nx2n matrix with A, B forming the first and the second 
n columns, respectively. This parametrization is obviously non-unique, since A, B 
can be replaced by CA, CB with any regular n x n matrix C. This defect can be 
corrected by choosing the matrices in the standard form [23124) . 

({/- J)#(0) + i(Z7 + I)#'(0) = 0, (2.3) 

where U is an n x n unitary matrix; the Hamiltonian corresponding to this condition 
will be labelled as Hjj . Elements of this family are labelled by n 2 real parameters 
which is, of course, the right number because all the Hjj are self-adjoint extensions 
of a common symmetric restriction with deficiency indices (n, n) [3]. 

Let us next recall a few examples of the boundary conditions (|2.3p . As men- 
tioned in the introduction, the requirement of continuity at the vertex selects a 
one-parameter subfamily corresponding to the so-called 5 coupling, 

n 

^(0)=V*(0) =:#>), j,k€n, £ ^(0) = c^(0) , (2.4) 

where a £ R and for brevity we have introduced the symbol n:={l,2,...,n}. We 
can add the case corresponding formally to a — oo, when the system decomposes 
into n halflines with Dirichlet endpoints, however, it is not interesting as long as 
we are concerned with nontrivial vertex couplings. In the particular case a = 
we speak about free boundary conditions since for the 5 function on line, n = 2, 
this corresponds to a free motion (sometimes the term Kirchhoff b.c, not very 
appropriate, is used). In terms of (|2.3|) the 5 coupling corresponds to the matrix 
U = — r — J — I, where J denotes the n x n matrix whose all entries equal one. 
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The 8' interaction on the line has two possible analogues for n > 2 25126] ■ One 
is a counterpart to (|2.4p called S' s coupling with the role of \& (0) , ^/'(0) interchanged, 

n 

$(o) = ^(o)=:V'(o), i.fceft, J>i(O) = 0V»'(O), (2.5) 

i=i 

where /3 e R U {+oo}. It corresponds to t/ = I — j^p-^i i n particular, the case 
/3 = oo refers to full Neumann decoupling. The other one, called 8' coupling, is 

n 8 
£>}(0) = 0, Vi(0)-^(0) = ^(^(0)-VU0)) J,fcen, (2.6) 

3=1 n 

with /? G K U {+00} which corresponds to U = + ^z^J- 

All the above examples have a common property, namely that the corresponding 
operators are invariant with respect to permutation of the edges, which is clear from 
the fact that matrices U are not changed by a simultaneous permutations of the rows 
and columns. The most general family of Hjj with this property is characterized 
by two parameters, U = al + bJ with \a\ — 1 and \a + nb\ — 1, cf. [20], the 
corresponding boundary conditions being 

( a -l)(^-(0)-Vfc(0)) + i(a+l)(^(0)-^(0)) =0, j,ken, 

(a - 1 + nb) 4>k (0) + 1(0 + 1 + nb) ^ ^(°) = ■ 
fc=i fe=i 



2.3. Approximation of 8 couplings 

Let us next recall briefly known results about approximations of vertex couplings 
starting from the 8 coupling. The idea is the same as for 8 interactions on the line. 

Let Us (a) :— n ^_ la J — I be the corresponding matrix of the condition (|2.3[) . 
Given a family of real- valued functions W — {Wj : j € n}, for simplicity assumed 
to be compactly supported, we define scaled potentials at graph edges by 

We,i := - e Wj(^) , 3 Eh. (2.8) 

Starting from the free boundary conditions and choosing the family (|2.8p we can 
approximate any nontrivial 8 coupling as the following result shows. 

Theorem 2.1. Suppose that Wj £ i 1 (0, 1) for j € n, then 

H Us{0 ) + W e —> H Us[a) as e -> 0+ (2.9) 

in the norm resolvent sense, where a :— 53?= i Jo 1 ^ / i( x ) ^ x • 

For proof see [TS] where a more general result of this type is derived, together with 
other extensions of the standard Sturm-Liouville theory to star graphs. 
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2.4. Approximation of singular permutation- invariant couplings 

Consider further permutation-invariant couplings with wave functions discontinu- 
ous at the vertex. Denote the operator Hjj corresponding to U — al + bJ with 
a, b satisfying the stated conditions as H a,b . The approximating family can be con- 
structed as follows: we start from the operator H Ui q := H Us ^ u ) and pass to H U)V 
obtained by adding a 6 interaction of strength v on each edge at a distance d from 
the centre. We will let the S's approach the centre scaling properly u,v. 

Theorem 2.2. Fix a pair of complex numbers a/-l and b ^ such that \a\ = 1 
and \a + nb\ = 1, and set 

n ( a — 1 + nb a — 1 \ 1 , ,. 1 a— 1 

u(d):=i—[ - + , v(d ):=---i -. 2.10 

w d 2 \a + 1 + nb a+lj w d a+1 y ' 

Suppose that a + l + nb ^ and a(a + nb) ^ 1, then the operators H u (d), v (d) converge 
to H a b in the norm resolvent topology as d — > 0+. Moreover, the claim remains 
true in the two excluded cases, provided we replace the above u(d) by —ndT x and 
£d~ u with R 9 ( / and v > 2, respectively. 

Proof can be found in [20], the particular case of S' s coupling (|2.5[) in which u(d) = 
—j3d~ 2 and v(d) = —d^ 1 was discussed in [19] . 

3. CS-type approximation of singular couplings 

After the preliminaries let us turn to our proper task, namely approximations of 
singular couplings a la Cheon and Shigehara, i.e. by means of additional 6 inter- 
actions, properly scaled, on edges of our star graph, without the requirement of 
permutation invariance. 

3.1. The class of approximable couplings 

The first question is how large is the class of operators Hjj which can can be treated 
in this way. We are going to answer it using the technique of [16], i.e. looking into 
convergence of the corresponding boundary conditions. 

Proposition 3.1. Let T be a star graph with n semi-infinite edges and T(d) be a 
graph obtained from T by adding a finite number of vertices at each edge. Consider 
a family {T(d) : d G R + } of such graphs with the properties that the number of the 
added vertices at each edge is independent of d and their distances from the centre 
are 0{d) as d — > 0+. Suppose that a family of functions ^ £ W 2 ' 2 (T \ ({c} U Vd)), 
where c is the centre ofT and Vd is the set of added vertices, satisfies the conditions 
\ 6 2,.1$ with d-dependent parameters, and that it converges to G W 2 " 2 {T \ {c}) 
which obeys the condition \2.1)) with some A, B satisfying the requirements i2.2\) . 
The family of the conditions \2.1}) which can be obtained in this way depends on 2n 
parameters ifn>2, and on three parameters for n = 2. 
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Proof. The 6 coupling in the centre of T is expressed by the condition (|2.4[) . 
Consider first 6 interactions on a halfline and look how the boundary values change 
when we pass between different sites. Suppose that at a point x the function and 
its derivative have the right limits, and that x + e is the site of a d interaction, then 
the Taylor expansion gives 

ip{x + e_) = 4>{x+) + eip'{x + ) + 0(e 2 ) , i//(x + e_) = tf(x+) + 0(e) , 

and the 5 interaction is according to 1|2.4|) described by 

ip(x + e + ) = ip(x + e_) =: ip(x + e) , tjj'(x + e+) - ip'(x + e_) = a(e)^(a; + e) , 

where a(e) is the coupling parameter. The may be e-dependent but we suppose 
such a dependence that the error terms can be neglected as e — > 0+; then we have 

</>(* + e) = 1>(x+) + etf(x + ) + 0(e 2 ) , 
V/(z + e+) = ^'(z+) + O(e) + a(e)(V(a!+) + ^'(*+) + <?(e 2 )) = 
= (1 + a(e)e)V'(ar + ) + a(e)V(a:+) + O(e) + «(e)G(e 2 ) , 

so that i/}(x + e) and ^>'(x + e+) depend on ^>(a;+) an d ip'(x+) linearly up to error 
terms. In case of a finite number of 6 interactions on a halfline one can show in a 
similar way recursively that the function value and the right limit of the derivative 
at the site of the last 5 depends, up to error terms, linearly on the function value 
and the right limit of the derivative for the first S interaction. 

Let us apply this conclusion to the edges of our star graph. We denote by dj 
the distance of the last 8 interaction on the j-th halfline family of edges in r(d); by 
assumption we have dj — 0(d). Then we have 

ifw^o = gf\d)m + hfm^) + rf\d) , 

If \dW j {d j+ ) = ~gf\d)m + hf\d)^(0) + ff\d) 

for some ff\ gf\ hf\ ff\ gf \ hf ) : K+ -» K. The functions ff ] and ff\d) 
are error terms and we suppose that they can be neglected in the limit. We are 
interested in the situation when the last relations can be inverted and i/)(0), i/)'j(0) 
can be expressed by means of ipj(dj) and ijj'Jdj,), 

m = f?\W>iVi) + gf ] (d)^(d J + ) + K(d) , j&h, (3.1) 
^ (0) = ff (dtyj (dj) + ~gf (dWj (d j+ ) + K{d) , j e n , (3.2) 

where we have introduced lZ(d) as the symbol for a generic remainder; we still 
assume that it can be neglected with respect to the other terms as d — > 0+. The 
equations (|3.1[) yield for j,k € n the conditions 

ff\d)i>j(dj)-$\d)^ h (d k ) + gf(d)^(d j+ )- j, ken 

(3.3) 
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and from (|3 - 2[) together with the second one of the conditions (|2.4[) we get 

n 

o^(0) = J] {} ( k\d)^{dk)+9t\dW k {d k+ )) + 11(d). 



k=l 



(3.4) 



We substitute for ^>(0) from (|3.ip and perform a repeated summation of (|3.4| over 
j. After an easy rearrangement we get 



E 



71 

(ajf(d) - nff V)) ^(d,-) + £ («5i X) (d) " ngf\d)) ^(d j+ ) = 71(d) 



3=1 



(3.5) 

Now we pass to the limit d — > 0+ in the equations ()3.3|) and (13. 5|) . Before that we 
multiply both sides by a power of d such that the right-hand side tends to zero 
as d — ► 0+, while at least one coefficient at the left-hand side remains nonzero, in 
other words, we use the assumed existence of the limit in which the error terms can 
be neglected w.r.t. the leading ones. The equation (|3.3|) acquires then the form 

(3.6) 



c^i(O) - c k M0) + Wj(0+) - t k ip' k (0+) = , j, k E 
while (|3.5p gives 



n 

E 

3 = 1 



7^(o) + E T ^K°+) = ° : 



(3.7) 



where Cj, tj, "/j, Tj are the appropriate limiting values of the functions involved. 
The obtained conditions can also be written in a matrix form, 



fa 


-C2 


•■ 


• ^ 




(h 


-t 2 


•• 


■ ° \ 


Cl 





-c 3 • ■ 


• 


*(0) + 


tl 





-is ■• 


■ 


Cl 





■■ 






h 





•■ 




\7i 


72 


73 •■ 


• In 1 






T-2 


r 3 . . 


■ T n ) 



*'(0) = 0. 



(3.8) 



It is clear already now - from the fact that the coefficients Cj , tj , jj , Tj, j G n are 
real- valued - that the achievable number of parameters cannot exceed An. 

So far we have not brought the self-adjointness into the game. To find the true 
number of parameters we pass from A, B to the unitary matrix U of standard 
boundary conditions (|2.3|) . This is achieved by multiplying the relation (|3.8p from 
the left by a regular matrix M such that U — I = MA and i(U + 1) = MB. This 
determines U since the last relations imply 

U=^M(A-iB), I = ~M(A + iB); 

notice that A + iB is regular because A and B are real and the matrix (vl|-B) has 
he full rank by assumption. Hence we have M — — 2(A + \B)^ 1 : which further gives 

U = -(A + iB)- 1 ■ (A-iB). 
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We shall apply the Gauss elimination method to get the chain of equivalences 
(-(A + iB)\{A - iB)) (I\ -(A + iB)- 1 ■ (A - iB)) ; 



the explicit form of A ± iB is obtained from (|3.8|) . We notice that the regularity 
of A + iB implies the following facts: (i) there is at most one j G n such that 
cj + itj = (and for such a j it holds that -fj + itj ^ 0) , (ii) there is at least one 



j G n such that jj + lTj ^= 



/-(d + rfi) 



C2 + 1*2 





The matrix (-(A + \B)\(A - iB)) equals to 



ci - i£i -(c 2 - it 2 ) 
c\ — i£i 



(ci+iti) ••• c n + it n ci-iti ••• —(c n —it n ) 

\ -(71 + lTl) -(72 + 1T 2 ) ' • • ~(7„ + 1T„) 71 - lTl 72 - 1T 2 ' ' ' 7„ - 1T„ / 

Suppose first that Cj + i*j ^ for all j G then by equivalent row manipulations 
we pass to the matrix (D\V), where 

/-(•n+m + fe+i)^®) 

C 2 + 1*2 

D = c 3 + 1*3 



/( 



V = 



ci 



V 



iti)S- 

2i cm 
S ci+ii 



2i cin— 7iti 
5 ci+iti 





O. ClTl— 71*1 

ci+iii 

■71*1 







C„ + l*n / 



-C 2 



1 C 2T2— 72*2 

1*2 



C2+l*2 

1 2l C2T2 —72*2 
h S C 2 +lt2 



_ 2l C»T»--y„tn 

Cn+l*n 
2l C n T n —^nt n 

~S 



»+i*r> 



2l C2T2 —72*2 
S C2+lt2 



1* + — 



2i c„t, 
S c„+it 



where we have denoted S = 53™= 1 Zl+vtl ' Since we used only equivalent manipula- 
tions, the diagonal matrix I? should have the same rank as A + i£>, hence it must be 
regular because none of its diagonal elements is zero. Consequently, we can divide 
each row of (£>|V0 by the corresponding diagobal element of D. This yields (I\U), 
where U is the sought unitary matrix and its diagonal and off-diagonal elements 
are given by 



u jk - 



McjTj - *j7j) _ c 3 - itj 

te + i^ELi^ c J + *i 

2»(c fc T fc - * fc 7fc) 

71 + 1T{ 



(3.9) 



if j^k. 



( Cj + i*j)(c fc + i* fe ) ELi 77+1*7 

The right-hand sides make sense due to the first of the conditions (|2.2[) and our 
assumptions about non- vanishing of all the expressions Cj + itj . 

So far we have not employed the second one of the requirements (|2.2p , namely 
the self-adjointness of the matrix AB* . This is equivalent to unitarity of U , however, 
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it is easier to check it in its original version. By a straightforward computation we 
find that the product AB* = AB T equals 



/ C\t\ + C 2 £ 2 C\t\ C\t\ 

c\t\ c\t\ + czh c\t\ 
c\t\ c\t\ c\t\ + Cit& 

C\t\ C\t\ C\t\ 

V7i*i _ 72*2 7i*i _ 73*3 7i*i - 74*4 



cm - c 2 r 2 \ 
cm - c 3 r 3 
cm - c 4 r 4 

cm - c n r„ 



Cl*l 
Cl*l 
Cl*l 

ci*i + c n t n 

71*1 - 7n*» 71 T l + 72T2 H 7n T n / 



hence A_B* is self-adjoint if and only if cm — CjTj = 7iii — T/tj holds for all 
j = 2, . . . , n, and therefore 



cm - 7i<i = c 2 t 2 - 7 2 i 2 = C3T3 - 73*3 = • • • = c „t„ - 7„i„ . 



(3.10) 



We denote the common value CjTj — jjtj as k and recall that we have denoted 
S = Yle=i c+iF 1 then the matrix {/ given by (|3.9p can be simplified, 



2lK 



(ci+iii^S ci+iti 

2lK 


( Cl +iti)(c 2 +it 2 )S 

2lK C 2 — 1^2 


(ci+iti)(c, 

2lK 




)S 


(c 2 +lt 2 )(ci+ltl)S 


(c2+lt2)^S C 2 +lt 2 


(c 2 +it 2 )(c, 




)S 


2lK 


2lK 


2lK 




c n —it n 


(c»+lt„)(ci+lti)S 


(c„+ii„)(c 2 +it 2 )S 


(c„+ii„)(c„+it„ 


)S 


c n +\t n 



Let us show that the matrix (|3.9|) can be parametrized by 2n real numbers. We 
rewrite the quantity S introduced above in the following way, 



/ / r 2 I +2 / < 



1 

''^2 r 2 , *2 ' 

fcl ^ + ^ 



and make first several observations: (i) regarding (|3.8j) as a system of linear equa- 
tions its solvability is not affected if the last one is multiplied by a nonzero number. 
At the same time, the value of k is directly proportional to jj, tj, and consequently, 
one can suppose without loss of generality that k = 1 (the case k = gives rise 
to the same situation as c\ + iti =0 which we shall discuss below), (ii) if k = 1 
the imaginary part of S is determined only by the values of Cj , tj, j £ n , (iii) and 
finally, one can also suppose without loss of generality that \c% +i*i| = 1, since in 
the opposite case we can divide all but the last of the equations in the system (|3.8[) 
by ci + iti| which is nonzero by assumption. 

With the above convention we can denote c% + itx =: & e and Re 5 —: p so that 



1 1 



E 

1=2 



1 

c# ~\~ tp 
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and U can be written explicitly as 



U 



/§ e" 2ie - 

2i 

(c 2 +it 2 )S 



2i 



(C2+1* 2 )S ^ 

2l C 2 — 1*2 



(c„+it„)5 
2i 



(c 2 +it 2 ) a S 



C2+l*2 



(c 2 +lt 2 )(c„+lt„)S 



v 

\ (o„+it„)S 



(c„ + lt„)(c 2 + lt 2 )S 



(o„+ii„)(c„+rf„)S 



C„+lt„ / 



(3.11) 



, *2,*3, ■ ■ -t n ,P- 



being dependent on 2n real parameters 9, C2, C3, . . . 

The above argument applies to any n > 2. In the case n = 2 the situation is 
somewhat different, because we have n 2 = 2n = 4 but (|3.11|) does not give the whole 
family of unitary 2x2 matrices; notice that the off-diagonal elements coincide. It 
is easy to show that the admissible U can be for n — 2 characterized by three real 
parameters. Indeed, writing U — (? b ) the unitarity requirement reads 



|6| 2 = 1, 



ld 2 = l, 



ab + be = . 



Knowing the modulus and phase of a, the modulus of b is determined so one has 
to choose its phase. Furthermore, since we assume 6^0 the element c is uniquely 
determined. Hence the matrix U of (|3 . 1 1 1) is described by three parameters which 
can be chosen, e.g., as the real parts of Ujj and the phase of U±2- 

Returning to the general case one can also write the conditions (|2.ip explicitly 
in terms of the parameters. A straightforward way is to put A = U — I, B = i(U + I) 
with U given by ([3. lip . To get a simpler expression one can pass from the system 
A 1 $(Q) + .B 1 i r '(0) = to an equivalent one multiplying it from the left by the matrix 



/ -& 9 c 2 + 1*2 
-& 9 c 3 + it 3 

-e ie 
V e w c 2 + it 2 c 3 + ii 3 



\ 




Cn + tin I 



this yields an explicit parametrization of the conditions (|2.ip with 



.4 



/ cos# 
cos 9 



cos 1 
Scos9 - 



-C2 







-C3 



Cl+ltl 



Sen 



B = 



sin 9 
sin 6 

sin 9 
Ssin9 1 1 



1 c„ 1 

c 2 +it 2 3 c 3 +it 3 



-ta 



us 



Stn -f- 7—7— St% + 7— 7— 

Cl+ltl z C2 + lt 2 ° C 3 +l«3 



C„+lt 



(3.12) 



Sir. 



C„+lt„ 
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and concludes the argument in the generic case when Cj + itj ^ for all j G n. 

It remains to deal with the case when the last mentioned requirement is violated; 
without loss of generality we may suppose that c\ + \t\ = 0. The matrix {—(A + 
iB)\(A - lB)) has then the form 



/ 







c 2 + it 2 





+ ITl) -(72 + IT2) 





C n -\- xtji 
-(in + 1T„) 







71 - in 



-(c 2 - it 2 ) 


72 - IT2 







( lt n ) 

7n - 1T n / 



Using the Gauss elimination scheme we arrive at (-D|U) with a diagonal D and 
upper-triangular V , and from here in the same way as above to (I\U) with 



C2 T 2 (-2*2 



u 



7l+ lr l 7l+ 1T l C 2 +lt2 

-^^XT 3 - 

C 2 +lt2 





V 







71+iTi C 3 +lt 3 



C3+l*3 







71+iTi c„+it„ 





c„+it„ 



Furthermore, it follows from the condition (|3.10| with c% = t\ = that 

C2T2 - 72^2 = C3T3 - 73^3 = • • • = C„T„ - 7„t„ = . 

hence all the off-diagonal elements in the above matrix U vanish which means that 
it is characterized by n real parameters, 

f/ = diag{e l9l ,...,e ie "} . 

It is easy to rewrite the boundary conditions in the form (|2.1[) and check that they 
correspond to the fully separated case, 



sin ■ 



^(o) + cos^-^(o) = o, jeh, 



2 TJ "~' ' 2 rj 

which is, of course, trivial for the viewpoint of quantum mechanics on T. 



(3.13) 
□ 



3.2. A concrete 2n-parameter approximation 

Knowing the maximum number of parameters in the boundary conditions (|2.1[) 
which can be achieved in this way, we are naturally lead to the idea of placing two 
5 interactions at each of the n halflines. In this section we are going to concretize 
this proposal. We will concentrate at the matrix (|3.11[) in the generic case leaving 
out the trivial situation (|3 . 1 3|) mentioned at the end of the previous proof. We will 
also leave out the case n = 2 which was discussed in the paper [27] , 

Let us specify the approximation arrangement. The <5's are placed as sketched 
on Fig. [TJ all dependent on a parameter d in terms of which the limit is performed: 

• there is a S coupling with parameter u(d) in the star centre 
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Approximations of vertex couplings in quantum graphs 13 

• on each halfline there is a S interaction with parameter Vj(d), where j is 
the halfline index, at a distance D(d) from the centre (it will turn out in 
the following that we may choose D(d) = d 3 ) 

• furthermore, each halfline supports another 5 interaction with parameter 
Wj(d) at the distance D(d) + d from the centre 




Fig. 1. Scheme of a 2n-parameter approximation 



For the sake of brevity we will not indicate the d-dependence of the parameters 
u, Vj, Wj and the distance D unless necessary. The boundary conditions which the 
functions ipi, . . . ,ip n on T have to satisfy are 

n 

V>i(0) = V> 2 (0) = ---=Vn(0)=:V(0), 5>K°+)=^(0) (3.14) 

1>i{D+) = i>j(D_) =: ipj (D) , $(£>+) - ^(D_) = v 3 ^{D) (3.15) 
i>j(D + d±) =: ipj(D + d) , 1/.(D + d+) - ^(D + d_) = w^j{D + d) (3.16) 

Further relations which will in the following serve to determine the parameter de- 
pendence on d are obtained from Taylor expansion of the respective wave functions, 

^(D) = ^(0) + D^(0 + ) + O(D 2 ) 7 $(£>_) = $(O+) + 0(£>), (3.17) 
i>j(D + d) = yjj(D) + #;•(£>+) + 0(d 2 ) , tf.(D + d_) = ip'j{D+) + 0{d) (3.18) 

for j G n. We need to find relations between the values tpi (D + d) , . . . , ip n (D + d) 
and tp[(D + d + ),..., tp' n (D + d + ). To this aim we express them first in terms of 
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0(0) and 0'-(O+). Using the relations (|3 . 1 5[) and (|3 . 1 7|) we get 

$(!>+) = $(0+) + O(D) + vM(0) + D^(0+) + 0(D 2 )) = 
= ^-0(0) + (1 + v j D)ip' j (0+) + O(D) + Vj O(D 2 ) . 

Substituting into the first one of the relations (|3.18[) and using (|3 . 15|) again we find 

ipjiD + d) = (1 + d^)0(O) + {D + d(l + v s D)) 0^(0+) 

+ 0{D 2 ) + dO{D) + d Vj O(D 2 ) + 0(d 2 ) . (3.19) 

The already obtained expression for ipj(D + ) together with the second one of the 
relations (|3.18j) give 

i/>'j(D + d_) = w 3 -V(0) + (1 + ^L>)Vj-(0+) + O(-D) + ^0(D 2 ) + 0(d) . 

Substituting from here and (|3. 19|) into the second one of the relations (|3.16|) we get 
after a simple rearrangement 

ip' j {D+d + ) = (vj + Wj (l + dvj)) 0(O)+(1 + VjD + Wj (D + d(l + Vj D))) 0^(0+) 
+ 0(D) + VjO(D 2 ) + O(d) + Wj (0(D 2 ) + dO(D) + d Vj O(D 2 ) + 0(d 2 )) . (3.20) 

Next we eliminate 0^(0+) (for simplicity we write 0^(0)) from the obtained rela- 
tions (f3~19|) and (|3^D| . multiplying them by 1 + VjD + Wj (D + d(l + VjD)) and 
D + d(l + VjD), respectively, and subtracting. In the resulting expression the coef- 
ficient at 0(0) equals one, 

(1 + Vj D + Wj {D + d(l + vjD))) ipj(D + d) 

= 0(0) + (D + d(l + VjD)) ip'j(D + d+) + Kj , (3.21) 

with the remainder term 

Kj := (1 + VjD + Wj (D + d(l + vjD))) (0(D 2 ) + dO(D) + dv 3 0{D 2 ) + 0(d 2 )) 

- (D + d(l + Vj D)) {(O(D) + v,0{D 2 ) + O(d) 

+ w 3 (0(D 2 ) + dO{D) + d Vj O(D 2 ) + 0{d 2 )) ) . 

So far the edge index has been kept fixed. Subtracting mutually the relations (|3.21[) 
for different values of j, k E N, we can eliminate 0(0), 

(1 + VjD + Wj {D + d(l + vjD))) ipj{D + d) 

- (1 + v k D + w k (D + d(l + v k D))) fe (D + d) 

= (D + d(l + VjD)) i>'j(D + d+) - (D + d(l + v k D)) 0^(D + d+) + Kj - K k 

(3.22) 

Returning to the relations (I3.19P and (|3.20j) we can eliminate from them 0(0) in a 
similar way as above arriving at the relation 

(1 + dvj)i/j'j(D + d+) - {Vj + Wj(l + dvj)) ipj(D + d) = 0^.(0) - Kj (3.23) 
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Approximations of vertex couplings in quantum graphs 15 

with the remainder term 

TZ 3 := ( Vj + wj{l + dv/)) (0(D 2 ) + dO(D) + d Vj O{D 2 ) + 0{d 2 ))- 
-(1+dvj) ((O(D) + v 3 0{D 2 ) + 0{d) + w s (0{D 2 ) + dO(D) + d Vj O(D 2 ) + 0{d 2 ))) 
Summing the above relations over j £ N and using (|3.14[) we get 

n n n 

^2{l+dv J )^ J (D+d + )~J2 ( V J + w 3^ + dv J^ ^A D + d ) = W>(0)-^ ilj . (3.24) 

3=1 i=l 3=1 

The right-hand side can rewritten using the continuity condition ()3 . 14[) in combi- 
nation with the relation (13.21|) . 

n n 

^(o) = I £ Vi(o) = I E ( o- + v i D + w i (° + d ^ + ^ D + rf ) 

71 3=1 n j=l 

- (D + d(l + VjD)) ip'^D + d+) + Tlj) . 
This allows us to cast (|3.24[) into a form which contains neither "0(0) nor ^-(0), 



71 

X { v i + w i( 1 + dv i) + - (! + V 3 D + w i ( D + d ( l + v 3 D ))i) V>j(-D + d) = 

n n 

= X (l + dvj + ^(D + d(l + Vj D))) ^(D + d+) + X - ■ (3-25) 

3 = 1 3 = 1 

The equations (|3.22|) and (|3.25|) are the sought relations between the function values 
and derivatives at the sites of the "outer" S's with ^(0) and ^ (0) eliminated. 

In the next step we are going to choose the dependences D — D(d), u = 
u(d), Vj = Vj(d) and Wj = Wj{d) for j £ h in such a way that the limit d — > 0+ 
will yield the (2n-parameter family of) boundary conditions (|2.1[) satisfying the 
requirement (|2.2[) . It appears that a suitable choice is the following one, 

D(d) := d 3 

1 a ■ 

1 + VjD — ctjd, i.e. Vj(d) 



rf 3 d 2 

1 + tOjd = [3jd , i.e. Wj(d) := — — + j3j 
u(d) := £ 



(3.26) 



(3.27) 



Indeed, in such a case the coefficients in p.22j) acquire the form 

(1 + VjD)(l + Wjd) + WjD = (otjPj - l)d 2 + Pjd 3 , 

D + d(l + VjD) = otjd 2 + d 3 

and a straightforward computation shows that the remainders are IZj = d 2 0{d), 
hence dividing (|3.22j) by d 2 we arrive at 

(otjpj - 1 + f3jd)ipj(d 3 + d) - {a k /3 k - 1 + (3 k d)^ k (d 3 + d) 

= (a, + d)^(d 3 + d+) - (a k + d)ij' k (d 3 + d + ) + 0{d) . 
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Taking the limit d — > + we have to realize that the condition ipj e VF 2,2 (R + ), j G h, 
requires that ipj(d) = o(d -1 / 2 ) holds at the halfline endpoint, hence we have 

{arfj - l)Vi(O) - (a k /3 k - l)^k(O) = «i^-(0) - «fcVfe(0) . j,k&n. (3.28) 



In a similar way we proceed with the equation (|3.25|) . We employ (|3.27|) . then a 
straightforward computation gives for the coefficients at ipj (D + d) and M (D + d + ) 
the following expressions 

u 

Vj + Wj(l + dv-j) + - (1 + VjD + wj (D + d(l + VjD))) 

= + <±( a .p. _ i)) i_ + (a. p. _ i + ^ 1 +0j , 

l + ( i, + ;(D + d(l + ^)) = (-l + ^)I + (a j + ^ + l, 

and the remainder terms IZj and — TZj are both <j _2 0(<i). We substitute from here 
to (|3.25j) . multiply the result by d 2 and pass to the limit d — > 0+; this yields 



3=1 



A - i)j ^(o) = E (- 1 + - a $(°) . ( 3 - 29 ) 

3=1 71 



The relations (|3.28p and (|3.29p are the sought boundary conditions. It remains 
to express them as (|2.ip and to find relations between the parameters contained in 
them to those of (|3.11[) . The matrix form of (|3.28p and (|3.29[) looks as follows, 



(aipx - 1 -(a 2 p2 - 1) 



otxPx - 1 
V 7i 





72 





-(a„f3„ - 1) 

In 



\ ( -a-y a 2 

*(0) - 







— ot\ • • • a, 

V n f 2 • • • f n / 



*'(o) = o. 



(3.30) 



where := -(aij/3j — l)—/3j and := 1 — —otj. We know that the corresponding 
matrix of (|2.3p is given by [/ = — (A + iB)^ 1 ■ (A — 1B), its matrix element being 



2i 



[a^j-l-iaj ^=1 (o!(ft _i)2 +Q 2 -w + i^ =1 (aiA _ 1) 2 +a 2j 

2i 

for j 7^ k. If the latter should correspond to (|3.1ip , it is sufficient to require 

\a1P1 - 1 -iai| = 1 (3.31) 



and 



U jk = 
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and to set 

EPi(ai/3i - 1) + ai 

|=1 WrDH^^ (3 - 32) 

ctjftj — 1 = Cj , —ctj = tj . (3.33) 

For a\ = the condition (|3.31[) is satificd trivially, while for a nonzero value it is 
equivalent to ct\ (a.\(ji\ + 1) — 2/3%) = 0, in other words we have to put 

In this way we have eliminated the parameter ai, and just 2n of them is left. 
The correspondence between the 2n-tuples fix, fa, fa, ... fa, oli, «3, . . . ,a n , u> and 
6*, c 2 , C3, . . . c n , i2, ^3, ■ ■ ■ > tn, P looks as follows: 

• fa < — > 0: they are related by |^ = e lS 

• otj,f3j < — ► Cj,tj, j 6 {2, . . . , rt}: see (|3.33p . 

• uj < — ► p: see (|3.32p . 

In what follows we will work with fa, fa, fa, ■ ■ ■ fa, 0.2, ot$, ... a n , u> , for simplicity 
we will use also a± remembering that it is determined by fa and the relation (|3.3ip . 



4. Norm-resolvent convergence 

The approximation worked out in the previous section was in the spirit of [16127) 
being expressed in terms of boundary conditions. One asks naturally what can 
be said about the relation between the corresponding operators. We denote the 
Hamiltonian with the coupling (|3.30p in centre of the star as H u ' a ^ , and H M ' v ' w (d) 
will be the approximating family constructed above, with a pair of 5 interactions 
added at each halfline. Our aim here is to demonstrate the following claim. 

Theorem 4.1. Let u, Vj, Wj, j € n, depend on d according to (|3.26p . i.e. 

Then H u ' v,w (d) converges to H^' '" in the norm-resolvent sense as d — > 0+. 

Proof. We have to compare the resolvents R Hu ^,nj^{k 2 ) and R H ^,sj(k 2 ) of the 
two operators for k 2 in the resolvent set. It is clearly sufficient to check the conver- 
gence in the Hilbert-Schmidt norm, 

\\R H ^ {d) {k 2 ) - R H ^ s j(k 2 )\\ 2 ^0+ as d — » 0+ , 

in other words, to show that the difference of the corresponding resolvent kernels 
denoted as G^' v ' w and G^'"' 13 , respectively, tends to zero in L 2 ((R + ) 2n ). Recall that 
these jernels, or Green functions, are in our case n x n matrix functions. 
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Let us construct first 0^' a '^ for the star-graph Hamiltonian referring to the 
condition (|2.ip in the centre. We begin with n independent halflincs with Dirichlet 
condition at its endpoints; Green's function for each of them is well-known to be 

K 

where x < := min{a; , y} , a;> := max{ir,z/}, and we put in = k assuming ReK > 0. 
The sought Green's function is then given by Krein's formula [H App. A], 



. , . ' ' L (( R+ )" 

where Rn{k 2 ) acts on each halfline as an integral operator with the kernel Q K and 
for 4>j(k 2 ) one can choose any elements of the deficiency subspaces of the largest 
common restriction; we will work with ^4>j(k 2 )(x)) = Sj m e~ KX . 

To find the coefficients Xji(k 2 ) we apply (|4~Tj) to an arbitrary * E L 2 (M+) 
and denote the components of the resulting vector as hj ; it yields 

/■ + OC n p + co 

hj(xj) = / g iK {x, yj )iP 3 { y] )A y] + V Xjiik 2 ) / e- Km ipi( yi )d yi ■ tT"*' . 
Jo l=1 Jo 

These functions have to satisfy the boundary conditions in the centre, 

n n 

£ A jm h m (0) + B jm h' m (0) = for all jeh. (4.2) 



= c KVm we find 

x m =0 



Using the explicit form of Gi K {x,y) and 

n I-+OQ 

MO) = VA m/ (fc 2 ) / e- Ky 'Myi)dyi (4.3) 

and 

r+oo n p-\-oo 

h' m (0)= e-^^ m (y m )dy m -Ky^\ ml {k 2 ) e^'^)^ ■ (4-4) 

Jo l=1 Jo 

Substituting from these relations into (|4.2j) we get a system of equations, 
" r+oc / n n \ 

£ / Z A?™ W* 2 ) + - K S B 3?nX m i(k 2 ) e~ Ky 'Myi)dyi = o , 

i=l ^° \m=l m=l / 

with j G ri. We require that the left-hand side vanishes for any ?/>i, i/>2, ■ ■ ■ , V'nj this 
yields the condition AA + £> — k£>A = 0. From here it is easy to find the coefficients 
Xji(k 2 ): we have (A — kB)A = —B, and therefore 

A iJ (* a ) = -[(A-«B)- 1 S] JI . 

Notice that the matrix A — is regular in view of the first conditions in (|2.2[) ; 
since A, B are real and Im k / 0, the requirement rank(A, i?) = n implies that we 
have also rank(A — kB) = n. 
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Let us now concentrate on the class of couplings for which we established in the 
previous section the boundary condition convergence. In this case A — kB equals 

( ai(/3i + K)-l -(a 2 (0 2 + K)-l) •■• ^ 

ai(/3i + K)-l •■• -{<Xn(j3 n +k) - 1 

V(/3 1 + K )(^a 1 -l)-^(/3 2 + K )(^ 2 -l)-^... (p n + K )^ an -l)-^J 

and a tedious by straightforward computation yields an explicit form of the matrix 
-(A - kB)- 1 B, namely 

[-(A- kB)~ 1 B] 
[-{A - kB)- 1 B] 



for j ^ i , 



In this way we get the Green function Qf^ a '^ ■ As we have mentioned above, it is an 
n x n matrix-valued function the (j, Z)-th element of which is given by 



1 1 



we use the convention that cc is from the j-th halfline and ?/ from the Z-th one. 

Next we will pass to resolvent construction for the approximating family of 
operators H u ' v ' w (d). As a starting point we consider n independent halfiines with 
Dirichlet endpoints; we know that the appropriate Green's function is Q 1K (x,y) = 
k^ 1 sinh nx < e~ KX> . The sought resolvent kernel will be then found in several steps. 
Each of them represents an application of Krein's formula. First we add the 5 
interaction with the parameter v at the distance d 3 from the endpoint, then another 
one with the parameter w at the distance d + d 3 , again from the endpoint. This is 
done on each halfline separately. In the final step we find Green's function for the 
star in which the Dirichlet ends are replaced by the 5 coupling with the parameter 
u. That will require, of course, to distinguish the halfiines by their indices. 

The first step is rather standard [193 an< i resulting Green function is 

S^(x,y) = g iK (x,y) - V G iK (y,d 3 ) g iK (x,d 3 ). (4.5) 

1 +v ■ £i K (u >« ) 

Adding another 5 interaction at the distance d from the previous one we seek the 
kernel in the form R v > w (k 2 ) = R v (k 2 ) + \{k 2 )((j)(k 2 ), -)(f>(k 2 ) where the first term 
is R v (k 2 ) := Q^ K and the deficiency-subspace element <fi(k 2 ) is chosen as 

0(k 2 )(x) :=g: K (x,d + d 3 ). 
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We apply this Ansatz to any feL 2 



and denote h := R v,w {k )ip. It is easy to 



check that Qf H (x,y) — G? K {x,y), hence we can write h explicitly as 

" + 0O 



h{x) 



gi{y,d + d s )i>{y)dyg^{x,d + d A ). 



By definition this function this function belongs to the domain of the operator with 
two 5 interactions, in particular, it has to satisfy the boundary conditions 



h(d + d 3 + ) = h(d + d 3 _) =: h(d + d 3 ) , 
h'(d + d 3 + ) - h'(d + d 3 _) =wh(d+ d 3 ) . 

Green's function continuity implies (14. 6j) . Furthermore, we have 



h'(x) 



(4.6) 
(4.7) 



dx 



i>(y)dy+\(k*) / g: K {y,d+d 3 )^{y)&y 



dx 



which allows us to express h'{d + d 3 + ) — h'{d + d 3 _). The first term obviously does 
not contribute to the difference, while the contribution of the second one simplifies 

= — 1 to the form 



£ dQ(x,y) 

m view of — a 

ox 



V+ 



9G(x,y) 
dx 



h'{d + d\) - h'{d + d 3 _) = -X(k 2 ) 



G: K (y,d + d 3 My)dy. 



To satisfy (14 .7p the coefficient A (A; 2 ) must obey the condition 



+00 



[X(k 2 ) + w + w\{k 2 )g^{d + d 3 ,d + d 3 )] g? K (y, d + d 3 )^{y) dy = 



for any ip € L 2 (R + ), where we have taken Green's function symmetry with respect 
to the argument interchange into account. Consequently, the square bracket has to 
vanish and we get the formula for the kernel with two S interactions, 



\ + wg^{d + d 3 ,d + d 3 ) 



g? K ( y ,d + d 3 ) gux.d + d 3 



(4.8) 

The remaining step will be more complicated because we are going to introduce 
a coupling between different halflines working this with matrix-valued functions. 
Our tool will be again Krein's formula which now takes the form 



n 

R H ^(k 2 ) = R H ,,4k 2 ) + ]T a^ 2 ) (<Mn-) 



J\Z=1 



L2((R+)~) 



<Mfc 2 ). 



where the functions (f>j(k 2 ) will be chosen as 



(<t> 3 (k 2 ){x)) m = 5 3 



dg^^(x,y) 
dy 



y=o 
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We apply this Ansatz to an arbitrary ^ = {tpx, . . . , i/'n} T and denote the elements 
of the resulting vector as hj , explicitly 



hj(x) 



+ OO 



1=1 



dx 



ipi (y) dy 



dg%' ws (x,v) 



x=0 



dy 



(4.9) 



where we have used Green's function symmetry and the fact that its complex conju- 
gation is equivalent to switching from k to R. As before the functions hi, h 2 , ■ ■ . , h n 
have to satisfy the boundary conditions expressing the S coupling in the star centre, 



h 1 (0)=h 2 (0) = --- = h n (0) =: h(0), 
h[(0) + h' 2 (0) + ■■■ + h' n (Q) = u ■ h{0) 



(4.10) 
(4.11) 



for any ipi, . . . , ip n S L 2 (M + ). Let us first express hj(0). The first term in the above 
expression does not contribute since QJi' '(0, y) = G? K (0,y) — 0. The second one 
contains the value of Green's function derivative which can be expressed using (|4. 



dg%' w *(x, y ) 



dy 



dg? K (x >y ) 



dy 



v=o 



w 



dg? K (y,d + d 3 ) 



i + wg? K (d + d 3 ,d+d 3 ) 



dy 



g^d + d 3 ). 



y=o 



The first term is obtained from (|4.5p together with the explicit form of the "free" 
kernel {?i K (x,y): we have 



dgu^v) 



dy 



in particular, ^g^p^ 



x=y=0 



i + v-g 1K (dz,d3) 

1. This further implies 



MM 3 ) 



dg%' w *(x,y) 



dy 



y=o 



l + v r g, K (d 3 ,d 3 ) 

(e- 



l + w r Q%(d + d*,d + <P) 



C - Kd ' g 1K ( x ,d 3 ) 

(d+d 3 ) 



i + vj-g^d^d 3 ) 



- Kd3 g M (d + d 3 , rf 3 )) • g% (x, d + d 3 ), (4.i2) 



in particular. 



dy 



x=y=0 



= 1 . Putting these results together, we can simplify 



the expression for the boundary values hj (0) as follows 

"UK* ! I 
1=1 Jo 



Mo) = EM fc2 ) 



My) d v- 



x=0 
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Now we can find what is required to fulfill the conditions (|4.10p . i.e. hj(0) = h m (0) 
for all j, m G n. This is true provided 



£(A j; (fc 2 )-A m z(fc 2 )) / 



Ox 



ipi(y)dy = 0, 



cc=Q 



holds for any n-tuple of functions tpi, . . . , ip n £ L 2 (K + ) which is possible if 

^jl(k 2 ) = A m ;(fc 2 ) for all j,m £n, I Eh, 

thus we can simplify notation writing A/ := \ji(k ) for a fixed / £ n. 

Values of the coefficients Ai, . . . , A n can be found from the remaining condi- 
tion (|4.1ip . To this aim we have to find explicit form of h'j (0) . It follows from the 
expression (|4.9[) for hj(x) that 



h'm = 



dx 



i=i 



x=0 

dG^ w '(x,y) 



Ox 



ipi(y) dy ■ 



d dQ^ Wi (x,y) 



x=0 



dx 



dy 



The boundary condition (|4.11|) then requires that the expression 



E 



dy 



u ■ X, 



Ox 



^i(y) dy 



vanishes for any ipx, . . . , ip n , and this in turn yields 

d (dQ^ Wi {x,y) 



X, 



E 

3=1 



dx 



dy 



x=0 



for all I £ n 



showing, in particular, that A; does not depend on Z, which means that all the 
coefficients Xji (fc 2 ) are the same and equal to the right-hand side of the last relation. 

Before specifying the expression in the square bracket let us write down the 
formula for the (j, l)-th component of the sought Green function: we have 



eg J' J (x, y ) 

dy 



y=0 



dQ^' w ' (x,y) 
dx 



u-T n — 

a l^m=l dx 



dy 



y=0 



(4.13) 



x=0 



The first derivative in the numerator was found in (|4.12[) and by Green's function 
symmetry the other one is given by the same expression, with y replaced by x. 



7, 2008 15:39 WSPC/INSTRUCTION FILE vertapp 



Approximations of vertex couplings in quantum graphs 23 



The same relation allows us to compute 
evaluate the quantity appearing in the square bracket above, 



l {x,y) 



, in particular, to 



dy 



y=0j 



i + w m -g^(d + d 3 ,d + d 3 ) 

I'm 



i + v m -g 1K {d\d?) 

^-K(d+d 3 ) 



■ e 



g 1K (d + d 3 ,d 3 ) 



(4.14) 



The relations (|4. 1 3[) and (|4.14[) together with (|4. 1 2[) and its mirror counterpart 
describe completely Green's function G^ v ' w of the approximating operators. 

After deriving explicit expressions for the resolvent we can pass to our proper 
goal which is to prove that the matrix-valued kernel G^l converges to G\ 
d — > 0+ which in terms of their components can be written as 



as 



lim 

d^0+ 



g U ,S0 



L 2 (R+xR+) 



= 0. 



Depending on the values x, y the difference G^ V j'i"( x i V)~ ^wli \ x > v) t a kes different 
forms. Notice that one can suppose without loss of generality that x < y, and 
therefore there are six different situations to inspect, namely 

• d + d 3 < x < y, 

• d<x<d + d 3 <y, 

• < x < d 3 , d + d 3 < y, 

• d 3 < x < y < d + d 3 , 

• < x < d 3 < y < d + d 3 , 

• < x < y < d 3 . 

To express the kernel difference we employ Taylor expansion of G^j^ix, y). Let us 
start with expressions which appear in the formulae repeatedly. The first one is 



l 

"3? 



d 2 



i + v m -g 1K (d3,d3) 

Using sinh(x) — x + 0(x 2 ) and e a 



sinh kc/ 3 c 



(*) 



1 + 0{x) we get 



sinh nd 3 



(Kd 3 + 0(d*))(l+0(d 3 )) 



K K 

and this in turn allows us to express (*) as follows, 
I 1 - a m d 

I* J 



= d 3 (l + 0(d 3 )), 



d 3 l + 



d? 



(d»(l + 0(d3))) 



1 



0(d) 



The next frequent expression is w m (l + w m ■ G^™ (d + d 3 ,d + d 3 )) . We employ 
relation (|4.5[) with v = v m and the expansion e x = 1 + x + 0(x 2 ) together with the 
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explicit form of Q 1K \ this yields after a straightforward computation 



g^(d+d 3 ,d+d 3 ) = d[ l-nd + 0(<f) =Q 1K (d+d 6 ,d+d d ) + 0(d 3 ), 



and therefore 



i + w m -g?,r(d + d 3 ,d + d 3 ) d 2 \ p m + K --L 



0(d) 



Now we can expand the first term in G^j'i (x,y). Using (|4.8|) for the parameters 
v — Vj, in = in j together with the previous result we get 



1 



sinh kx e Ky 



1 



d 2 yj3 m + K -J- 
sinh c~ K!/ 1 



0(d) 



sinh «;(d + d 3 ) e KJ/ sinh «(d + <i 3 ) e 



/3m + « " 



(l + 0(d)) 



As for the second term in f (|4.13|) ). we first expand the derivative in the denominator 
using Q 1K (d + d 3 , d 3 ) = d 3 (l + 0(d)) and (gHU). A direct computation yields 



dx \ 



dy 

and therefore 

d / egsr^{x,y) 







y=o) 


-H 


x=0 



0m ~l~ K 



0(d) 



E 



d.?; 



dy 







y=oj 






x=0/ 



En 
m—1 a r 



0m+K 



•0(d) 



,(/3 m +K)-l 



Next we expand the derivatives which appear in the numerator using the relation 
g?™ [x, d + d 3 ) = d(l + 0{d)) e~ KX ; it gives 



dg? r Wm {x,y) 



dy 



i + v m -g 1K (d 3 ,d 3 ) 

Wrn 

i + Wm -g?^(d + d 3 ,d + d 3 ) 
1 



e~ Kd ' g iK (x,d 3 )- 



d 2 \ (J3 m + k) - 1 



0(d) e 



and the analogous expression for 



with x replaced by y. This de- 



termines the behaviour of the second term at the right-hand side of ((HUH])) as 
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d — > 0+, and for the full kernel G^ V ji"(x, y) we consequently have 



/sinh/tx e~ Ky 1 + 0(d) 



V K ft ■ + « - ^ 



+ 



+ 0(d) e-^e-^ . 



On the other hand, for x < y we have 



+ 



sinh ki e Ky 



1 1 



e e 



hence the Green function difference satisfies 



Cr(^)"CA^)=0(<i)e- ra e-«» as d^ 0+ . 



The same estimate is obviously valid also for d < y < x, hence there is a constant 
K independent of d, x and y such that 



< Kde- KX e~ Ky 



(4.15) 



holds for all d < 1, x > d + ci 3 and y > d + d 3 . Now we are in position to estimate 
the Hilbert-Schmidt norm of the resolvent difference for the operators H^' '" and 
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H u ' v ' w (d) which can be written explicitly as follows, 



\R H u.^ (d) (k 2 )- R H »,a,ft\\ 2 - 

3,1=1 



OO f + OO 



o Jo 



dxdy 



3,1=1 



OC f+OO 



d+d 3 J d+d 3 



r>u,v,w ( \ s-ruj,a.8 / \ 



d+d 3 r +oc 



d 3 J d+d 3 
+oo r d+d 3 



+ 



d+d 3 Jd 3 

d 3 r+oo j 



J d+d 3 
+oo r d 3 



'( x >y) - G?K°jf( x >y) 



d+d 3 Jo 
d+d 3 r d+d 



GikJi fay) - g^j fay) 



dxdy 
2 

dxdy 

2 

dxdy 
dxdy 

2 

dxdy 



d 3 Jd 3 
d 3 r d+d 3 



*lK,jl 



fay) -y*^°ji ( x ,y) 



dxdy 



o Jd 3 

d+d 3 r d 3 



d 3 JO 



y^A ( x >y)-G^.u ( x >y) 



y X nM ( x >y) -SiK.n fay) 



dxdy 



dxdy 



o Jo 



y^M fa v) - GikA ( x ' v) dxd v 



The inequality (|4.15[) makes it possible to estimate the first one of the integrals, 



+oo p+oo 



d+d 3 J d+d 3 



y,nji ( x >y)-yiK,ji fay) 



2j2 



dxdy < K 2 d 



e -2(ReK)x dx 



d+d 3 



= K 2 d 2 



3 -2(Re K)(d+d J ) 



2Re k 



< 



K 2 



2Rck 



and it is obvious from this inequality that for d — > + the integral tends to zero for 
any j, I £ n. In a similar way one can estimate each of the remaining eight integrals: 
using Taylor expansions of Q^ji we § e ^ a bound for the integrand which shows 
that the integral vanishes as d — > 0+. Since the argument repeats the procedure 
described above, we skip the details. Putting all this together, we conclude that 

Jim \\R H u,e,vs {d) (k 2 ) -R Ha , 3 j\\ 2 = 0, 

and therefore the resolvent difference tends to zero in Hilbert-Schmidt norm as 
d — > 0+ which is what we set up to demonstrate. □ 
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5. Approximations with added edges 

We have seen that a CS-type scheme can produce a 2rt-parameter family of (self- 
adjoint) couplings out of the whole set depending on n 2 real numbers. To get a 
wider class we have to add to the star graph T not only vertices but edges as well. 

5.1. Admissible couplings 

The first question naturally is how many parameters can be achieved in this way. 
An upper bound on this number is given by the following statement. 

Proposition 5.1. Let T be a star graph with n semi-infinite edges and denote by 
{f (d) : d E M+} a family of graphs obtained from T by adding finite edges connect- 
ing pairwise the halflines; their number may be arbitrary finite but independent of d. 
Suppose that T(d) supports only S couplings and 5 interactions, their number again 
independent of d, and that the distances between all their sites are 0(d) as d — > 0+. 
Suppose that a family of functions ^ £ W 2 " 2 (T \ ({c} U Vd)), where c is the centre 
ofT, and Vd is the set of the vertices added on the halflines, satisfies the conditions 
with d-dependent parameters, and that it converges to "J £ W 2 " 2 (T \ {c}) 
which obeys the condition V2.1\) with some A, B satisfying the reguirements \2.2]) . 
The family of the conditions \2.1\) which can be obtained in this way has real-valued 
coefficients, A, B G M n,n , depending thus on at most ("J 1 ) parameters. 

Proof. The S coupling in the centre of r(d) , identified with centre of T, is expressed 
by the conditions (|2.4[) . For any j € n we denote by dj the coordinate of the most 
distant point on the j-th halfline which supports either a 5 interaction or a 5 
coupling at the endpoint of an added edge. We arrange the function values at these 
points into the n-tuple ^(d), and similarly ^'(d + ) is the n-tuple of right derivatives. 
Let us stress that this a symbolic notation; the elements are tpj(dj) and i// (ci ) +), 
respectively. 

As in the proof of Proposition 13.11 we can use ()2.4|) to express these quantities 
through the common value tp(0) and the right derivatives 4 , '(0+) at the origin 

Mi(d)*(d) = i/>(0) ■ m 2 {d) + M 3 (rf)*'(0) + 11(d) , 

JVi(d)*'(d+) = 4>(0) ■ n 2 {d) + JV 3 (d)*'(0) +n{d) 

for some Mi, M 3 , N t , N 3 : R+ -> R"*™, m 2 ,n 2 : R+ -> M™ and error terms K,K : 
M. + — > 1" supposed to be negligible as d — > 0+; we may assume that 1Z, 1Z = o(l). 
The above system can be also written in a matrix form, 

'MM -m 2 (d) -M 3 (d) \ 

Ni(d) -712(d) -N 3 (d) J +> = 

- Q 11 - 1 '\n6)) 

To find an approximation in the described sense one has to find a relation between 
*&(d) and \t'(d+) eliminating ^>(0), ^'(0). Since the former are determined by the 
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latter we may suppose that the matrices M\{d) a N\{d) are regular; the elimination 
then leads to a system 

A(d)*(d) + B(d)*'(d+) - K(d) , 

where the matrices A(d), B(d) are real for all d G M + and the right-hand side 
consists of an error term TZ : K + — > M n . We multiply the last equation by a power 
of <i such that the right-hand side is o(l) as d — > 0+ while the left-hand one has 
a nontrivial limit. It is clear that we can get in this way the condition (|2.3p with 
real- valued coefficients, A,BG M. n ' n . □ 

5.2. A concrete approximation arrangement 

The above discussion leaves open the question how such an approximation can be 
constructed to cover the mentioned (™^" 1 )-parameter family. Our aim here is to 
demonstrate a specific way to do that. We consider the coupling (12. 3|) with real 
A, B, and for simplicity we restrict our attention only to the generic case assuming 
that B is regular so that the boundary conditions acquire the form 

*'(0) = -B _1 A*(0) 

with a symmetric matrix —B~ 1 A. We can also write them as 

V(0) = (D + S)*(p), (5.1) 

where the real matrix D is diagonal while S is real symmetric with a vanishing 
diagonal; it is clear that D and S depend on n and (J) real parameters, respectively. 

To construct approximation of the corresponding operator H A ' B we have find 
suitable family of graphs T(d). The decomposition of the matrix in (|5.1jl into the 
diagonal and off-diagonal part inspires the following scheme: 

• the centre of T supports a 5 coupling with the parameter u{d) the depen- 
dence of which on d will be specified below 

• at each edge of T we place a 5 coupling at the distance d from the centre; 
the corresponding parameter Vj(d), to be again specified, will be related to 
the diagonal element Djj of the matrix D 

• the pairs of edges whose indices j, k correspond to nonzero elements of 
the matrix S we join by an additional edge, whose endpoints are the S 
coupling sites mentioned above, and in the middle of this edge we place the 
S interaction with a parameter w^j ^ (d) related to the value of Sjk 

The metric on T and T(d) is intrinsic, nevertheless, it is useful to think of it as 
of induced by embedding of the graphs into a Euclidean space. Without loss of 
generality we may consider the original star T as a planar graph and to construct 
as embedded into R 3 . In such a case, of course, we have to make sure that the 
added edges do not intersect. This can be achieved in the way sketched in Fig. O 
A possible way is to employ the bijection b from the family of two-element subsets 
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Fig. 2. Approximating graph: a star amended by connections of the edges, with a <5 coupling in 
the centre, one 5 coupling at each edge and one <5 interaction at each (broken) connection segment 



of {1, 2, . . . 7i } to the set {1, 2, . . . "C" -1 ) }. The edge connecting the j-th and fc-th 
halflinc is formed by two segments connected in a V-shape. Its endpoints are at 
the j-th and fc-th halfline, both at the distance d from the centre. The tip of this 
1^-graph is placed on the halflinc starting from the centre of T in the perpendicular 
direction to its plane - see Fig. [3]- at the distance bjk -d 2 , so that the length of the 
connecting V-graph is dyT + (bjkd) 2 . 



As before we denote by tpj the wave function on the j-th halfline assuming that 
all the coordinates have zero in the centre of T. Furthermore, we denote by ip^ 

and <pjf' k ^ the wave function on the line segment part of the connection between 
the j-th and fc-th halfline which is attached by one of its endpoints to the j-th and 
fc-th halfline, respectively; notice that the order of the upper indices is irrelevant. 
Such a connecting link is regarded as a star with two edges of the same length. For 




Fig. 3. The connecting edge between the j-th and fc-th halflinc 
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the sake of brevity we introduce also the set Nj defined as 

N 3 ; = {k e n : S jk + 0} ; 

its cardinality #Nj tells us how many nonzero elements are in the j-th row of the 
matrix S, in other words, how many V-shaped connecting edges sprout from the 
point Xj — d on the j-th halfline. 

Next we will write down the boundary conditions describing the involved S 
couplings; for simplicity we will not indicate the dependence of the parameters 
u, Vj, W{j.fc} on the distance d. The S coupling in the centre of T means 

n 

^(0)=^ 2 (0) = --- = ^ n (0)=:^(0), ^VK°+)=^(0), (5.2) 

3=1 

the S interaction at the "tip" of the broken edge connecting the j-th and k-th 
halfline between the vertices added at the distance d from the centre (of course, for 
j,k 6 n such that Sj k ^ only) is expressed through the conditions 

v { / M (o) = =■■ ^ fe} (o) , (^ k} y J (o + )+(^ k} y k (o+) = W{m ^ k} (o) , 

(5.3) 

and finally, the S coupling at the mentioned added vertices added requires 

^ (d + ) = ^ (d_ ) = ^- k] (<yi + (b jk d) 2 ) =: ^ (d) , j e h , k e 

^(d+) - ^(d_) - £ fceJVj (^ J '' fc} )'(dVl + - «iifc(d) , jsn. (5 - 4) 

Further relations which will help us to find the parameter dependence on d come 
from Taylor expansion, 

^(d) = yj J (0) + d^(0) + O)(d 2 ), ^.(d_)=^-(O + ) + 0(d), J En, (5.5) 

vf' fe} (dy/l + (b ]k d) 2 ) = ^(0) + d^l + {b ]k d) 2 {^< k] )'(0 + ) + 0(d 2 ) , 

V , . _ x (5-6) 

{^■■ k] y Uyji + (b jk dy_) = (^/■ k} )'(o + ) + o(d) , j, k £ a, 

where we have used the fact that y/l + {bjkd) 2 = 1 + 0(d 2 ). Now we employ the 
first of the relations (|5.6|) together with the continuity (15. 4p , which yields 



d^l + {b jk d) 2 (^ m )'(0 + ) = ^(d) - v {3M (0) + 0{d 2 ) . (5.7) 

The same relation holds with j replaced by k, summing them together and using 
the second of the relations (|5.3[) we get 



2 + dy/l + {b jk df w {m j <p^ k H0) = ^-(d) + Vfe(d) + 0(d 2 ) . 

We express (p^ ,k '(0) from here and substitute into (|5.7p obtaining 

^ /i , ^ ^2 / {i,*)\'/n \ , / ,n ^(d) + ^ fc (d) + 0(d 2 ) 2 

d\fl + (b jk d) 2 (ipy ') (0+) = Vj(d) - ,,/■,,,, no + C(d). (5.8) 



February 7, 2008 15:39 WSPC/INSTRUCTION FILE vertapp 



Approximations of vertex couplings in quantum graphs 31 

The relations (|5.5p and (|5.2[) give 

d^(0+)-^(d)-^(0) + O(d 2 ), (5.9) 
and summing this over j E ft we arrive at the identity 

n n 
3=1 3=1 

The right-hand side of it can be rewritten using (|5.2p . This makes it possible to 
express ip(0); substituting it into (|5.9|) we get 

#ko + ) = - ELi t^r (d2) + °^ ^ 

Next we use consecutively the second relations of (|5.4j) . (|5.5[) and (|5.6[) to infer 

^(d + ) =V^(d) + (^>)' fc (d^l+^d^ J + $(d_) 

=«^i(d) + £ (^^(0+) + $(°+) + 0(d) . 

keNj 

Substituting into the last relation from (|5.8|) and (|5.10j) we get 



V'j(d) + Md) 



y/1 + (b jk d) 2 2 + dyjl + (b jk d) 2 ■ w {:j , k} 



where we have also employed the fact that O(d) [1 + [bj^d) 2 ]^ 1 / 2 = 0(d) holds as 
d — > 0+ for all j 7^ fc, j, k £ n. 

Now we can finally ask about the parameter dependence on d. Since the last 
relation is supposed to yield in the limit d — > 0+ the j'-th row of the matrix condition 
(|5.ip . it would be sufficient to have the following requirements satisfied: 

lim [v j +- 1 \ V ; + 1 - V 1 I I = D 

<^°+V V^Vi + (W 42 + vi + (M) 2 % t} jj 

(5.11) 

for all j G n, 



3 



lim i • - 1 1 = (5.12) 

<i-o + d yiTTvdF 2 + dy/i + (b Jk d) 2 w {3 , k} 



for all j 7^ fc, j, k £ h, and finally 

1 



d(n + du) ^ ^ ^ ^ 
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as d — + 0+. To fulfil (|5 . 1 2[) one can choose 

which makes sense because Sjk ^ by assumption, since then the limit equals 

i™+ 1 + O(cP) 2d + (1 + Q(d 2 )) (-^ - 2d) ~ ' 

With the choice (|5 . 14[) taken into account the condition (|5.1ip will be satished 
provided Vj + g(#iVj + 1) - Efeew, S jk = Dj, i-e. 

Vj .(d):=^.-#^-±l- £ (5.15) 

Finally, the last requirement will be satisfied, e.g., if the expression equals d which 
is true if 

Summarizing the argument we conclude that choosing the parameters in the de- 
scribed approximation according to (|5.14jl — (|5.16j) we get in the limit the generic 
boundary conditions (|5.1[) . We conjecture that such an approximation would again 
converge in the norm-resolvent topology. 
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